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Abstract 

Using the tadpole improved clover Wilson quark action on small, coarse and anisotropic 
lattices, KN scattering length in the 1 = 1 channel is calculated within quenched 
approximation. The results are extrapolated towards the chiral and physical kaon 
mass region. Finite volume and finite lattice spacing errors are also analyzed and a 
result in the infinite volume and continuum limit is obtained which is compatible 
with the experiment and the results from Chiral Perturbation Theory. 
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1 Introduction 



Anisotropic lattices and improved actions have been extensively used in re- 
cent lattice QCD calculations. In our previous works, the following tadpole 
improved gluonic action on anisotropic lattices has been utilized: 
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where Pqi and Pij represents the usual temporal and spatial plaquette vari- 
able, respectively. Rij and Rio designates the 2x1 spatial and temporal Wilson 
loops, where, in order to eliminate the spurious states, we have restricted the 
coupling of fields in the temporal direction to be within one lattice spacing. 
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The parameter Us-, which is taken to be the fourth root of the average spa- 
tial plaquette value, implements the tadpole improvement. With the tadpole 
improvement factor, the renormalization of the anisotropy (or aspect ratio) 
^ = tts/at would be small. In this work, we will not differentiate the bare 
aspect ratio and the renormalized one in this work. Using this action, glue- 
ball and hadron spectra have been studied within quenched approximation 
[1,2,3,4,5,6,7]. In our previous works, configurations generated from this im- 
proved action have also been utilized to calculate the tttt scattering lengths in 
the 1 — 2 channel [8]. In this paper, we extend our study to the kaon nucleon 
scattering lengths within quenched approximation. Unlike the case in pion 
pion scattering where the chiral expansion is quite reliable, chiral predictions 
in KN scattering suffer from large corrections [9,10,11,12]. On the experimen- 
tal side, Martin used dispersion relation fits to determine the KN scattering 
parameters [13]. Therefore, a lattice QCD calculation will offer an important 
and independent check on these results. 

Lattice calculations of hadron scattering lengths have been performed by var- 
ious authors using conventional Wilson fermions on symmetric lattices with- 
out the improvement [14,15,16]. In Ref. [15], a systematic study of the hadron 
scattering was performed, including pion-pion, kaon-nucleon, pion-nucleon and 
nucleon-nucleon scattering. However, this study was done only at some spe- 
cific parameters and the results were not extrapolated towards the continuum 
and chiral limit. In later works, chiral and continuum limit was studied sys- 
tematically in the case of pion-pion scattering, especially in the 1 — 2 channel 
[16,8,17,18] It is known that, using the symmetric lattices and Wilson action, 
large lattices have to be simulated which require substantial amount of com- 
puting resources. In Ref. [8], we have shown that, pion-pion scattering lengths 
can be obtained using relatively small lattices with the help of tadpole im- 
proved Wilson quarks on anisotropic lattices. It is also possible to perform the 
chiral, infinite volume and continuum limits extrapolation. The final extrapo- 
lated result on the pion-pion scattering length in the 1 — 2 channel is in good 
agreement with Chiral perturbation result, the experimental result and other 
lattice results using conventional Wilson fermions [17]. In this paper, similar 
techniques are applied and the kaon-nucleon scattering lengths in the I — 1 
channel is studied. 

The fermion action used in this study is the tadpole improved clover Wilson 
action on anisotropic lattices [19,20,8]: 
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where various coefficients in the fermion matrix A4 are given by: 
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Among the parameters which appear in the fermion matrix, csw is the coef- 
ficient of the clover term and u is the so-called bare velocity of light, which 
has to be tuned non-perturbatively using the single pion dispersion relations 
[20]. Tuning the clover coefficients non-perturbatively is a difficult procedure. 
In this work, the tadpole improved tree- level value, namely csw = 1, is used. 

In the fermion matrix (2), the bare quark mass dependence is singled out 
into the parameter a and the matrix A remains unchanged if the bare quark 
mass is varied. This shifted structure of the matrix A4 can be utilized to 
solve for quark propagators at various values of valance quark mass ttlq (or 
equivalently k) at the cost of solving only the lightest valance quark mass 
value K, — /^max, using the so-called Multi-mass Minimal Residual (M^R for 
short) algorithm [21,22,23]. This is particularly advantageous in a quenched 
calculation since one needs the results at various quark mass values to perform 
the chiral extrapolation. 

This paper is organized as follows. In Section 2, the method to calculate KN 
scattering lengths is briefly reviewed. The KN correlation functions used in 
our calculation are defined. In Section 3, some simulation details are described. 
In Section 4, our results of the KN scattering lengths obtained on lattices of 
various sizes and lattice spacings are extrapolated towards the chiral limit. 
Physical strange quark mass is also determined from the X-input. Finite size 
effects are then studied and the infinite volume limit are taken. Finally, our 
lattice results are extrapolated towards the continuum limit. Comparisons 
with the experimental value and vahics from chiral perturbation theory at 
various orders are also discussed. In Section 5, we conclude with some general 
remarks. 



2 Formulation to extract the scattering length 

In order to calculate the elastic scattering lengths for hadron-hadron scattering 
on the lattice, or the scattering phase shifts in general, one uses Liischer's 
formula which relates the exact energy level of two hadron states in a finite 
box to the clastic scattering phase shift in the continuum. In the case of KN 
scattering at zero relative three momentum, this formula amounts to a relation 
between the exact energy E^^\f of the KN system in a finite box of size L with 
isospin I, and the corresponding scattering length in the continuum. This 
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formula reads [24] : 
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where ci = —2.837297, C2 = 6.375183 are numerical constants and jiKN — 
rriKiTiN / {rriK + mjy) is the reduced mass of the KN system. In this paper, the 
KN scattering length q'q^ in the 1=1 channel will be studied. 

To measure the hadron mass values rriK, "tun and to extract the energy shift 
SE^jpj^, we construct the correlation functions from the corresponding operators 
in the 1 = 1 channel. In this work, the hadron states that need to be measured 
include pion, kaon, rho, nuclcon and KN states. The local operators which 
create a single pion with appropriate isospin values are: 



7r+(x,t) = -ci(x,t)75M(x,t) , TT (x,t) = ti(x,t)75(i(x,t) , 

7r°(x, ^) = ^ [""(x, t)j5u{yi, t) - J(x, i)75(i(x, t)] , (5) 

where w(x, t), d(x, t), ■u(x, t) and (i(x, t) are basic local quark fields correspond- 
ing to u and d quarks, respectively. The operators which correspond to zero 
momentum single pions are: 

<W = X^E^'^(X'^) ' (6) 



where the flavor of pions a = +, — , and is the three volume of the lattice. 
Zero momentum one pion correlation function can then be formed: 

aw =< <(tK(o) > . (7) 



Similarly, one can construct the kaon correlation functions from kaon opera- 
tors. In the 1=1 channel, only and proton operators are needed: 



X+(x,i) = -s(x,i)75ii(x,i) , 

p(x, t) = Sabc kr(x, t) (C~^75)4(X, t)] Uc(x, t) 



(8) 



The proton operator given above is itself a Dirac spinor. The zero momentum 
proton-proton correlation function is obtained by projecting out the positive 
parity part of the zero momentum spinor operator and taking the trace in 
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Dirac space: 

CN{t) - E 0) p{y, t)) . (9) 



Prom the large t behavior of the correlation functions, the pion mass m^, the 
kaon mass uik and the nucleon mass mjv are obtained. 

Two hadron operators arc used in the construction of two hadron correlation 
functions. For example, the KN operator in the 1=1 channel is given by: 

0'£Mt) = K^{t)po{t + 1) , 0'£j^{t) = K^{t)Po{t + 1) , (10) 



where KqII) and po{t) are the zero momentum kaon and proton operators 
respectively. From these operators, two hadron correlation functions can be 
constructed. The KN correlation function in the 1=1 channel reads: 

C'K=Nit)=<Ol,=^{0)(^)o^^=^{t)> . (11) 



Numerically, it is more advantageous to construct the ratio of the correlation 
functions defined above: 



n'=\t) = cirM/{CK{t)Cj,{t)) . (12) 



This ratio thus exhibits the following asymptotic behavior for large t: 

7^^=^(^)*>^^-<-* , (13) 



with SEj^'j^ = Ej^'j^—mK—'mM is the energy shift in this channel which directly 
enters Liischer's formula (4). This formula is utilized for small SE^j^j^^t values 
(large enough L) in the calculation where the signal to noise ratio is good. In 
this case, one could equally use the linear fitting function: 

7^^=l(^)^»^»^l-5£;gj^^ , (i4) 



to determine the energy shift 6E\^\^. Similar procedure has also been used in 
lattice studies of pion-pion scattering lengths by various groups. 

KN correlation function, or equivalently, the ratio TZ^^^it) constructed above 
can be transformed into products of quark propagators using Wick's theorem. 
These contributions are represented by diagrams as shown in Fig. 1, where the 
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(a) 



(b) 



Fig. 1. The diagrams that contribute to KN scattering in the 1 = 1 channel. The 
soUd hnes represent u-quark propagators while dashed lines are d-quark propagators. 
The dashed-dotted lines designates the strange quark propagator. Note that diagram 
(b) in fact contains two different contractions. 

solid, the dashed and the dashed-dotted lines represent the u, d and s quark 
propagators, respectively. ^ The KN correlation function in the 7 = channel 
is, however, more complicated which involves vacuum diagrams that require 
to compute the quark propagators for wall sources placed at every time-slice, 
a procedure which is more time-consuming than in the 1 = 1 channel. 



3 Simulation details 

Simulations are performed on a PC based cluster with 30 nodes. Configura- 
tions are generated using the pure gauge action (1) for 6^40, 8'^40 and lO'^SO 
lattices with the gauge coupling P = 1.7, 1.9, 2.2 and 2.4. The spatial lattice 
spacing is roughly between 0.2fm and 0.4fm while the physical size of the 
lattice ranges from 1.2fm to 4.0fm. For each set of parameters, several hun- 
dred decorrelated gauge field configurations arc used to measure the fermionic 
quantities. Statistical errors are all analyzed using the usual jack-knife method. 
Basic information of these configurations is summarized in Table 1. 

Quark propagators are measured using the Multi-mass Minimal Residue al- 
gorithm for several different values of bare quark mass. Periodic boundary 
condition is applied to all three spatial directions while in the temporal direc- 
tion, Dirichlet boundary condition is utilized. In this calculation, it is advan- 
tageous to use the wall sources which enhance the signal [14,15,16,8]. Values 
of the maximum hopping parameter Kmax, which corresponds to the lowest 
valance quark mass, are also listed in Table 1. Typically, a few hundred Min- 
imal Residual iterations are needed to obtain the solution vector for a given 
source. On small lattices, in particular those with a low value of /3, the hop- 
ping parameter has to be kept relatively far away from its critical value in 

^ In this work, the u and d quark are assumed to be degenerate in mass. 
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Table 1 

Simulation parameters for lattices studied in this work. Input aspect ratio parameter 
^ is fixed to be 5 for all lattices being studied. The approximate spatial lattice 
spacing in physical units as obtained from [2,7] is also indicated. Also listed are 

the maximum value of the hopping parameter Kmax- In the last three columns, we 
listed the physical up/down quark hopping parameters, the physical strange quark 
hopping parameters and the results for the factor F = Oq^ m?p/ fiKN at physical 
quark mass values, respectively. In the rows labelled as "Extrapolated", we have 
also listed the values of F after the infinite volume extrapolation. 



Lattice 


No. confs 






r^s 


F = a^Q^my^KN 


P = 2.4, ui = 0.4090, as ~ 0.22fm, i/=0.93 


6^40 
8^40 

10^50 


288 
192 

204 


0.0610 
0.0610 
0.0605 


0.06179(2) 
0.061791(9) 
0.061636(9) 


0.06120(6) 

0.06113(4) 
0.06113(5) 


-0.7(15) 

-1.6(6) 
-3.3(9) 


Infinite 










-3.3(11) 


P = 2.2, ui = 0.378, as ~ 0.27fm, z/=0.95 


6^40 
8^40* 
10=^50 


256 
224 
184 


0.0605 
0.0600 
0.0605 


0.06116(1) 
0.06113(1) 
0.061230(7) 


0.06020(2) 
0.06041(5) 
0.06024(2) 


-1.9(3) 
-3.1(9) 
-3.3(13) 


Infinite 










-3.9(11) 


p = 1.9, u'l = 0.3281, as =i 0.35fm, iy=0.93 


6^40* 
8^40* 
10^50* 


304 
160 
252 


0.0600 
0.0595 
0.0595 


0.06087(1) 
0.06073(1) 
0.060722(6) 


0.05945(2) 
0.05969(5) 
0.05964(2) 


-2.3(7) 
-2.8(46) 
-4.0(44) 


Infinite 










-4.1(46) 


P = 1.7, = 0.295, ~ 0.39fm, z/=0.90 


6340* 

8^40* 
10^50* 


288 
240 
250 


0.0590 
0.0585 
0.0580 


0.06027(1) 
0.06019(1) 
0.06007(2) 


0.05895(3) 
0.05906(4) 
0.05914(9) 


-2.9(25) 
-3.1(58) 
-3.4(62) 


Infinite 










-4.6(64) 



order to avoid the appearance of exceptional configurations. The parameter 
z/, also known as the bare velocity of light, that enters the fermion matrix (2) 
is determined non-perturbatively using the single pion dispersion relations as 
described in Ref . [20] and Ref . [8] . The optimal value of v which corresponds 
to each parameter set of the simulation is also tabulated in Table 1 
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Single pion, kaon, rho and nucleon mass values are obtained from the plateau 
of their corresponding effective mass plots. The fitting interval is automatically 
chosen by minimal per degree of freedom. Due to the usage of finer lattice 
spacing in the temporal direction, good plateau behavior was observed in these 
effective mass plots. Therefore, contaminations from excited states should be 
negligible. These mass values will be utilized in the chiral extrapolation. 



0.9 




T 



Fig. 2. The ratio R^^\t) as a function of t obtained from the KN correlation func- 
tions for 8 different values of the bare quark mass. The lattice size is 8^40 and the 
gauge coupling /3 = 2.4. The optimal bare velocity of light is taken to be = 0.93. 
The straight lines are linear fits according to Eq. (14), from which the energy shifts 
5E = E^jpj^ — rriK — itin are extracted for all possible values of Kg)- 

KN correlation functions and the ratio 7^*^^^ (t) are constructed from products 
of suitable quark propagators according to Wick's theorem. Both the Direct 
and the Cross contributions are included. For the ratio TZ^^^t), we obtain 
clear signal for all our data sets. Since a kaon consists of one up (down) and 
one strange quark whose mass is different from the up and down quarks, all 
physical quantities involving the kaon will depend on two quark masses. In 
terms of the hopping parameters, the KN correlation functions depends on a 
pair of hopping parameters: {ku,Ks)- In Fig. 2 we show the ratio TZ^^^t) for 
different pairs of k values as a function of the temporal separation t, together 
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with the corresponding hnear fit (14) for /3 = 2.4 on 8^40 lattices. The starting 
and ending positions of the straight Hnes indicate the appropriate fitting range 
which minimizes the per degree of freedom. 



After obtaining the energy shifts SEj^j^, these vahics are substituted into 
Liischer's formula to solve for the scattering length Oq^^ for all possible hopping 
parameter pairs: k,,), which corresponds to the up (down) and the strange 
bare quark mass values respectively. This is done for lattices of all sizes being 
simulated and for all values of p. Prom these results, attempts are made to 
perform an extrapolation towards the chiral, infinite volume and continuum 
limit. 



4 Extrapolation towards chiral, infinite volume and continuum lim- 



The chiral extrapolations of physical quantities involving both the up (down) 
and the strange quarks consists of two steps. In the first step, the bare strange 
quark mass, or equivalcntly the corresponding hopping parameter Ks, is kept 
fixed while the hopping parameter of the up/down quark, Ku, is brought to 
their physical value k^'*^). The precise value of 

can be obtained by 

inspecting the chiral behavior of the pseudo-scalar (pion) and the vector me- 
son (rho) mass values. In the second step, one fixes the up/down quark hop- 
ping parameter at its physical value obtained in the first step and extrapo- 
late/interpolate in the strange quark mass. The physical strange quark hop- 
ping parameter k^'*^) is determined by demanding the mass of the kaon being 
exactly at its physical value (493McV). In practice, this could be accomplished 
by, for example, inspecting the kaon to rho mass ratio. This is also known as 
the i^-input. After these two steps, we would be able to obtain physical quan- 
tities at the physical up/down and strange quark mass values. In Table 1, we 
have hsted the physical hopping parameters k^'^^^ and k^'^^^ for the up/down 
quarks and the strange quark. Note that one only needs the pion, kaon and 
rho correlation functions to determine the physical hopping parameters of the 
quarks. 

We now come to the chiral extrapolation of the scattering length. In the chiral 
limit, the KN scattering length in the I — 1 channel is given by the current 
algebra result: 



its 



(1) 




(15) 



where a^Q^ is the KN scattering length in the I — 1 channel and fx 



~ 113MeV 
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is the kaon decay constant. Chiral Perturbation Theory to one-loop order gives 
an expression which includes the next-to-leading order contributions [12]. Un- 
like the chiral perturbation theory results for the mesons, the chiral expansion 
for baryons also receives corrections which are suppressed by the meson to 
baryon mass ratios. Since the mass of the kaon is not light, one typically ob- 
tains large corrections from the next chiral order compared with the previous 
one. In the case of the KN scattering length aQ^"*, the next order correction 
is 114% with the opposite sign and the next-to-next order is about 94% with 
the same sign of the tree-level. This is a signal that one really has no control 
over the expansion. 

To perform the chiral extrapolation of the scattering length, it is more conve- 
nient to use the quantity F = OQ^m^p/ hkn-, which in the chiral limit reads: 

F = ^ ^ = --r^ ^ -3-686 , 16 



where the final numerical value is obtained by substituting in the experimental 
values for rUp ~ 770MeV and fx ~ 113MeV. The factor F can be calculated 
on the lattice with good precision without the lattice calculation of meson 
decay constants. The error of the factor F obtained on the lattice will mainly 
come from the error of the scattering length aQ'\ or equivalently, the energy 
shift SEjpj^ and nowhere else. Since we have calculated the factor F for several 
different values of valance quark mass, we could make a chiral extrapolation 
and extract the corresponding results in the chiral limit. 

First, we extrapolate the factor F in the up/down quark mass values. In 
this step, we have adopted a simple linear extrapolation, the fitting range of 
the extrapolation is self-adjusted by the program to yield a minimal per 
degree of freedom. Second, the extrapolated values for the factor F are then 
extrapolate/interpolate in the strange quark mass. Our hopping parameters 
are chosen such that the physical strange quark mass is either within the 
range of our choice or very close to avoid a long extrapolation in the strange 
quark mass. In Fig. 3, we have shown the two step chiral extrapolation for 
one of our simulation points. In the lower panel of Fig. 3, we have plotted 
the results for the factor F for all possible pairs of quark mass values. They 
are arranged according to different up/down quark hopping parameter. The 
straight lines show the extrapolation in (the first step extrapolation). The 
physical up/down quark hopping parameter are read from the pion and rho 
correlation function fits which are tabulated in Tab. 1. The resulting factor 
F after the first extrapolation are plotted in the upper panel of Fig. 3 versus 
the strange quark mass parameter 1/(2ks). The straight line is the linear 
extrapolation/interpolation (second extrapolation) and the final result is also 
depicted as a solid square at the physical strange quark mass. Current algebra 
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Fig. 3. Chiral extrapolation for the quantity F 



rni/jiKN for our simulation 



results at /3 = 2.4 on 8^40 lattices. In the lower panel, the factor F as defined in 
Eq. (16) are plotted versus the pion mass squared, for fixed strange quark mass. The 
straight lines represent the corresponding linear fit in m\ for the data. In the upper 
panel of the plot, the extrapolated quantity F is plotted (open squares) as a function 
of strange quark mass parameter 1/(2ks). The straight line represents the linear 
interpolation/extrapolation to the physical strange quark mass, where the result is 
depicted as a solid square. As a comparison, the corresponding semi-experimental 
result [13] for this quantity is shown as a filled circle. Chiral perturbation theory 
results are also shown as a filled triangle (tree), a filled square (next) and a filled 
hexagon (next-to-next) , respectively. 



result and the result from chiral perturbation theory are also shown as a filled 
triangle (tree level), a filled square (next order) and a filled hexagon (next- 
to-next), respectively. The result from Ref. [13] is shown as a filled circle. To 
avoid crowdedness, these symbols are shifted to the left. The fitting quafity 
for the factor F are reasonable. The quafity for the chiral extrapolation of 
other simulation points are similar. As is seen, the linear fit gives a reasonable 
modelling of the data. The divergent contributions and other higher order 
corrections from quenched chiral perturbation theory seem to be numerically 
small for the lattices being simulated. 
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After the chiral extrapolation, we now turn to study the finite vohime effects 
of the simulation. According to formula (4), the quantity F obtained from 
finite lattices differs from its infinite volume value by corrections of the form 
However, the situation might be different in a quenched calculation, 
which happened in the case of pion-pion scattering [15,25]. In the case of KN 
scattering length, a similar calculation using quenched chiral perturbation 
theory is still not available. Therefore, we will be using a linear extrapolation 
in as suggested by Liischer's formula. In Table. 1, we have listed the 

infinite volume extrapolation for our simulation points at P — 2.4, 2.2, 1.9 
and 1.7. 
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Fig. 4. Continuum extrapolation for the quantity a'^^rrip/fJ,} 



obtained from our 



simulation results at /? = 2.4, 2.2, 1.9 and 1.7. The straight Une represents the linear 
extrapolation in a^/ro. The extrapolated results are also shown, together with the 
experimental result from Rcf. [13] indicated by the filled circle. For comparison, 
the results from Chiral perturbation theory are also shown as the points at = 
respectively. 

Finally, we can make an extrapolation towards the continuum limit by elim- 
inating the finite lattice spacing errors. Since we have used the tadpole im- 
proved clover Wilson action, all physical quantities differ from their continuum 
counterparts by terms that are proportional to Ug- The physical value of ag 
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for each value of [3 can be found from Ref. [2,7], which is also included in 
Table 1. This extrapolation is shown in Fig. 4 where the results from the chi- 
ral and infinite volume extrapolation discussed above are indicated as data 
points in the plot for all 4 values of P that have been simulated. The straight 
line shows the extrapolation towards the Ug — limit and the extrapolated 
results are also shown as solid squares together with the experimental result 
from Ref. [13] which is shown as the filled circle. For comparison, the results 
from chiral perturbation theory ar various orders arc also shown as the filled 
triangle, the filled square and the filled hexagon at = 0, respectively. It is 
seen from the figure that our extrapolated result is negative, with a large error 
bar. Therefore, we can only say that, from our preliminary study we obtain a 
result which is compatible with the experimental result and results from chiral 
perturbation theory. The error of our calculation can be reduced by further 
studies with more statistics and on lattices with smaller lattice spacing in the 
future. To summarize, we obtain from the linear extrapolation the following 
result for the quantity F = a^o^m^/ hkn = (—1-6 ± 5.4). If we substitute in 
the physical values, we obtain the KN scattering length in the 1 = 1 channel: 
a'o'^ = (— 0.17±0.59)fm, which is to be compared with the experimental result 
of (-0.33)fm and the chiral result of (-0.398 + 0.457 - 0.389)fm at the first, 
second and third chiral order, respectively. 



5 Conclusions 



In this paper, we have calculated kaon nucleon scattering lengths in isospin I = 
1 channel using quenched lattice QCD. It is shown that such a calculation is 
feasible using relatively small, coarse and anisotropic lattices. The calculation 
is done using the tadpole improved clover Wilson action on anisotropic lattices. 
Simulations are performed on lattices with various sizes, ranging from 1.2fm to 
about 4fm and with four different values of lattice spacing. Quark propagators 
are measured with different valence quark mass values. These enable us to 
explore the finite volume and the finite lattice spacing errors in a systematic 
fashion. The infinite volume extrapolation is made. The lattice result for the 
scattering length is extrapolated towards the chiral and continuum limit where 
a result consistent with the experiment and the Chiral Perturbation Theory is 
found. We believe that, using the method described in this exploratory study, 
more reliable results on KN scattering lengths could be obtained with more 
statistics and with simulations on finer lattices. 
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